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A qubit can be used as a sensitive spectrum analyzer of its environment. Here we show how the 
problem of spectral analysis of noise induced by a strongly coupled environment can be solved for 
discrete spectra. Our analytical model shows non-linear signal dependence on noise power, as well 
as possible frequency mixing, both are inherent to quantum evolution. This model enabled us to 
use a single trapped ion as a sensitive probe for strong, non-Gaussian, discrete magnetic field noise. 
To overcome ambiguities arising from the non-linear character of strong noise, we develop a three 
step noise characterization scheme: peak identification, magnitude identification and fine-tuning. 
Finally, we compare experimentally equidistant versus Uhrig pulse schemes for spectral analysis. 
The method is readily available to any quantum probe which can be coherently manipulated. 



The ability of a quantum system to withstand noise is 
characterized by its decoherence rate; the rate at which 
superpositions deteriorate. Hahn's discovery pQ of the 
echo technique showed that reducing decoherence can be 
achieved by external modulation, e.g in the case of spins, 
performing a single spin flip during the experiment. Since 
then, the idea of using external modulation to prolong 
coherence, known as Dynamic Decoupling [2JE], has been 
well developed to include many pulses [SJ [7J and different 
modulation schemes [5lll2j. 

These techniques rely on the condition that a quantum 
system, modulated at frequency /, is most influenced by 
the noise power spectral density at /. For a two-level 
system experiencing phase noise which is either Gaussian 
or weak (perturbative) , this condition takes the following 
integral overlap form p~3l415j : 

1 f°° 
R(T) = -- lnpu = / dfS(f) \F T (f)\ 2 (1) 

where p is the system density matrix, R(T) the deco- 
herence rate, S(f) the noise power spectral density and 
Fr(f) is the Fourier transform of the applied modulation 
for an experiment of length T. By generating modulation 
schemes which have a small overlap with the noise spec- 
trum, R diminishes. Therefore, different dynamical de- 
coupling schemes are optimal for different noise spectra. 
For a detailed account, see [H] and references therein. 

Can one use decoherence as a measurement tool? For 
that matter, Eq. ([I]) can be considered in the inverse 
manner. Instead of prolonging coherence by separating 
the modulation and noise spectra one can focus Fr(f) 
about a given spectral component at say, fo, thereby ex- 
tracting S(fo). This idea of a qubit-based spectrum ana- 
lyzer has been suggested in the context of different qubit 
technologies [17- 20 and has been recently analyzed for 
different types of noise spectra [23] ■ 

Experimental realizations of spectral analysis through 
spin decoherence spectroscopy were performed with dif- 
ferent technologies; for example, in trapped ions [22], cold 
atomic ensembles [23, 24J, nitrogen- vacancy centers in di- 
amonds [25] , super-conducting flux qubits [22] and NMR 



experiments in molecules [27] . The decoherence spectrum 
can be used to fix parameters in a known noise model, or 
reconstruct it through the inversion of the decoherence- 
spectrum relation [211 [27] . 

The operation of such qubit spectrum analyzers re- 
lies on the validity of Eq. ([I]), i.e. a linear response in 
the spectrum as in the case of weak or Gaussian noise. 
Otherwise, the qubit evolution can be significantly non- 
linear in Hamiltonian terms. In which case, the relation 
between measured decoherence and noise is very hard to 
calculate. 

It turns out, as will be shown in this paper, that spec- 
tral analysis in the strong noise limit can be significantly 
simplified for discrete spectra. This is reminiscent of 
the use of simple frequency analysis tools which enabled 
Babylonian astronomers to accurately predict the timings 
of lunar and solar eclipses [2E] and 19 th century scholars 
to provide tide predictions for various coasts and har- 
bors [29] . This is despite the fact that nonlinear evolution 
is present in planet and ocean dynamics as well. 

The merit of using a strongly coupled qubit can be un- 
derstood in terms of the trade-off between signal sensitiv- 
ity and spectral resolution using the Cramer-Rao bound. 
For an experiment time T and noise amplitude N the 
weak assumption is equivalent to a small noise index 
■q = NT <C 2-7T. If coherence is estimated by a quantum 
projective measurement, the optimal amplitude signal- 
to-(projection)noise ratio of a continuous spectrum is ob- 
tained when rj — > 0, consistent with the weak assumption. 
For large N, however, this will impose a short experi- 
ment time thereby limiting the spectral resolution which 
scales as 1/T. Better spectral resolution requires depar- 
ture from the weak limit. 

For the case of discrete spectra, the Cramer-Rao bound 
implies that both the amplitude signal-to-noise ratio and 
the spectral resolution are optimal for rj > 2ir. More- 
over, the spectral resolution attains an enhancement fac- 
tor and scales as l/(T^/fj) due to the non- linear response 
of the coherence with respect to noise amplitude. Spec- 
tral analysis of discrete spectra should therefore benefit 
from operating in the strong noise limit. 



2 



In this work we describe a simple analytical model ex- 
tending Eq. ([I]) to the non-perturbative, non-Gaussian 
discrete case. We further show how one can use this the- 
ory to identify the noise spectral components and mea- 
sure their magnitude in typical noise scenarios. We apply 
our spectral analysis scheme, using a single-trapped ion, 
to analyze the discrete spectrum of magnetic field noise 
in our lab. 

We focus on a two level quantum probe described by 
\ip(t)) = a |t) + ||) and governed by a Hamilto- 
nian H — h(N(t)& z + Q(t)a x )/2 where N(t) is classical 
dephasing noise and f2(t) is the spectrum analyzer mod- 
ulation. We assume no spin relaxation processes. Our 
purpose is to use the modulation tt(t) to quantify the 
noise N(t). 

For a probe initialized to |V>o) = (|t}+ |4-))/v / 2 the 
superposition relative phase at time T is |30j : 

pT poo 

<t>{T) = / dtN(t)F(t) = / dfN(f)F T (f) (2) 

JO J-oo 

where F(t) = cos(/* dt'Cl(t')), N(f),F T (f) are the re- 
spective Fourier transforms, the latter calculated on a 
truncated experiment window of length T . 

Phase coherence is obtained by averaging over noise 
realizations, A = (e 1 ^). The assumption that noise is 
either weak or Gaussian translates into (e l< ^) = e - ^ ^ 2 . 
Here, the decoherence rate R is proportional to (</> 2 ), the 
variance of the superposition phase imposed by noise. 
Combined with Eq. pi), relation (JlJ follows. 

To calculate the phase coherence without assum- 
ing that noise is Gaussian or perturbative, we as- 
sume discreteness: N(t) — J2k=i\Nk\cos(2nfkt + 
afe). For a single noise component, according to ([2]), 
<t> = \N Q F T (f Q )\ cos(a) so (e*) = (27T)" 1 C dae^ = 
Jo(|-^o-fr(/o)|) where Jo is the zeroth Bessel function 
of the first kind. For an ideal sinusoidal modulation at 
fo, A — Jo(|iVfc | T). The weak limit is valid only when 
Jo can be well approximated to second order in its ar- 
gument. Moreover, when the noise index rj = |A^ |T 
crosses z$ s» 2.4, the first zero of Jo (a;), coherence be- 
comes negative, i.e. the phase superposition partially 
refocuses close to tt. Notice that in general rj = c \Nk 1 T 
where c is a numerical constant depending on the mod- 
ulation shape. For the square wave modulation c = 2/n 
and for Uhrig modulation c m 0.42. Such single Bessel 
behavior due to a single mechanical resonance of a can- 
tilever coupled to a nitrogen-vacancy (NV) center was 
recently observed [3T] . 

In the case of more than one noise component, the 
coherence behavior takes a product form over all noise 
components. Assuming a k € [0, 2ir] are uniform mutu- 
ally independent variables, 

A(T)^(e^)=l[jo(\N k F T (f k )\) (3) 

k 



This equation coincides with Eq. ([T]) by Taylor expand- 
ing the Bessel functions to second order and recalling 
that R oc — In A Equation ([3| is the main tool of our 
noise spectral estimation method. 

The strong noise limit also reveals frequency mixing if 
we allow correlations in the a k -s. Whenever an integer 
combination of the noise frequencies is nulled J2k^kfk = 
0, additional Bessel product terms affect the coherence, 

d 

hi, hd k — 1 

Zh k f k =0 
S/ifc even 

(4) 

where the h k are integers and Jh k the corresponding 
Bessel functions of first kind. The dominant summand 
corresponding to hi = . . . = hd = coincides with Eq. 
Q. By focusing the modulation at a single frequency, 
as in our experiment, all the higher Bessel terms can be 
neglected, and information on the phase relation between 
different spectral components is lost. One will be able to 
retrieve it via the cos(S/ifeafe) term by using a multi-tonal 
modulation. 

Our system is comprised of the two spin states of 
the electronic ground level of a single 88 SV + ion, |f 
) =|5*i/2,J = lA Mj = 1/2) and ||) =\5s 1/2 ,J = 
1/2, Mj = —1/2). This Zeeman sub-manifold is first or- 
der sensitive to external magnetic fields. The dominant 
noise we measured was magnetic field fluctuations B(t) 
due to power line harmonics rendering a discrete noise 
spectrum, N(t) — g[iBB{t)/h, where g is the Lande g- 
factor, hb the Bohr magneton and H the Planck constant 
divided by 2tt (see Fig. [1^) . We performed spin rotations 
by pulsing a resonant rf magnetic field. Rotation angles 
and rotation axes were controlled by tuning the pulse 
duration and the rf field phase, (f> r f, respectivley. State 
initialization and measurement were performed by opti- 
cal pumping and state-selective fluorescence correspond- 
ingly 1331133]. 

To measure the phase coherence we performed a 
Ramsey-type experiment as shown in Fig. [lja. A mod- 
ulation f2(t) of length T is sandwiched between two tt/2 
pulses, differing by a relative phase <f> r f. We then mea- 
sured the probability of the ion to be in the |f) state P+ 
as a function of <j) r f . A fit to P-j- = | — ^ cos <f> r f yields 
an experimental estimate of the phase coherence A. Ex- 
amples of such fringes are shown in Fig. [lj: and [TJi. In 
all cases, £l(t) was a train of it pulses at different times 
and possibly different rotation axes. 

A first distinctive characteristic of non-perturbativity 
is the negative values of the coherence (noise index 
rj > Zq), shown in [TJ;. Here we fixed the modulation 
frequency at f m od = lOOPz while increasing n, the num- 
ber of equidistant pulses. As seen, the fringe contrast 
with n = 19 (shown in Fig. [TJi) is inverted with respect 
to n = 1 (shown in Fig. [lj:). A fit to Eq. ^ is shown 
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FIG. 1. (a) 8S Sy + ground state manifold comprising the 
Zeeman sensitive quantum probe, (b) Typical experimental 
sequence with n modulation pulses sandwiched between two 
7r/2 pulses with a relative <j> r f phase. (c,d) Probability of 
finding the probe in the |t) state at the end of a sequence 
vs. the rf relative phase. Red line is a best fit to P-\ = 
1/2 — A/2 cos((f) r f), A is the coherence. Plot lc was measured 
with n — 1 pulses and Id with n = 19 equidistant pulses with 
inter pulse distance of 5 ms. (e) Coherence vs. number of 
equidistant pulses n. The first and last points correspond to 
the fringes in lc and Id. Red line is a single parameter best-fit 
to Eq. pi. Inversion of the fringe contrast in d correspond to 
negative coherence, (f) Coherence vs. modulation frequenct, 
with n = 11 equidistant pulses. Although the spectrum shows 
five spectral features, it corresponds to a single, highly non- 
perturbative noise component. Red line is a single parameter 
best-fit to Eq. {3}. 



by the red line, assuming a single spectral component at 
100 Hz, with Nq as a single fit parameter and results in 
B 100 Hz = 3.0(2) M G. 

A second mark of non-perturbativity is that multi- 
ple spectral features can arise from a single noise com- 
ponents, as shown in Fig. [If. The number of pulses 
is fixed at n = 11 and the modulation frequency is 
scanned around / = 100 Hz. The spectrum shows a 
"power broadened" spectral feature around 100 Hz with 
five coherence minima. Unlike the perturbative case 
these do not correspond to five different spectral com- 
ponents but rather to a broadened response to a mag- 
netic field monotone. Again Eq. ^ with Nq as a sin- 
gle fit parameter shown by the red line is used, yielding 
-Sioo Hz = 15.3(3)/iG. This noise amplitude corresponds 
to a noise index of ?y = 10.3(2), well in the strong noise 
regime. The noise amplitudes extracted from the data 
shown [IJb and[lj are very different as these data sets were 
taken at different times. 

To practically estimate a multi-tone discrete spectrum 
we first identify the frequencies fk of its components. 
In any modulation scheme, the peak of the modulation 
Fx(f) increases linearly with the total experiment time 



T while improving spectral resolution. To identify the 
different noise components we therefore modulated the 
probe at different frequencies. For each modulation fre- 
quency the number of pulses was increased until the dif- 
ferent noise components emerged. Examples are shown 
in figures [2k and [2]d. Here, the measured coherence is 
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FIG. 2. Spectral peak identification, (a) Two scans taken 
several months apart. Each scan is a colormap of coherence 
vs. the number of pulses n and the modulation frequency. 
The upper scan shows clear features at 100 Hz and 150 Hz. 
The latter nearly vanishes in the lower scan and a new 200 Hz 
component appears, (b) Magnitude extraction of a single 
noise component. The number of pulses n is varied while the 
modulation frequency is kept fixed at 100 Hz until a zero- 
crossing is reached. From this zero crossing the field magni- 
tude can be estimated. Red line is a fit to Eq. |3| with a 
single fit parameter (c) Same as b for the 200 Hz component. 

shown vs. the (equidistant) modulation frequency and 
the number of ir pulses. As the number of pulses was 
increased clear spectral features emerged. The two data 
sets were measured four months apart with a different 
magnetic environment; the spectral response at 150 Hz 
which is clear in Fig. [2k almost vanishes in [2}) where a 
new 200 Hz component appeared. 

Once the component frequencies {fk} have been deter- 
mined, the multiplicative structure of equation ^ is used 
to determine their magnitudes Nk- Whenever the coher- 
ence A(T) crosses zero, with high certainty, only one of 
the Bessel functions in the product is nulled. If the mod- 
ulation is centered about fk, increasing the experiment 
time T until the first zero crossing occurs implies that 
the corresponding Bessel has been nulled and provides 
an estimate for Nk- We used this method to extract the 
magnitudes of the 100 Hz and 200 Hz components iden- 
tified in Fig. [2]d. We focused an equidistant modulation 
at / = lQQHz (200Hz) while increasing the number of 
pulses, n. Coherence vs. n is shown in Figure [2J; (|2ji) . 
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From the zero crossing we estimate a noise amplitude of 
2.6(2) /iG (5.0(2) /zG). This is in reasonable agreement 
with a best-fit to Eq. ([3| shown by the red line, yielding 
2.9(l)/xG (4.9(3) fj,G). 

The last stage of spectral characterization is fine tuning 
of the estimated noise magnitudes with a full fit proce- 
dure, using the previously estimated field magnitudes as 
a starting point. Such a fit to Eq. ([3| is shown in figure 
[3] with five fit parameters B 50 h z = 2.0(1)/jG, B WO h z = 
15.4(4)/iG, B 150Hz = 4.2(3) M G, B 2Q0Hz = 6.3(3)/xG and 
a slowly varying field, {g^BB s i ow /h)f s i ow = 66(2)Hz 2 . 
The non-perturbative nature of the spectrum is quan- 
tified by the corresponding noise indices: 7750 = 2.7(1), 
Vwo = 10.4(3), 77 15 o = 1-9(1), 77200 = 2.1(1). 
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FIG. 3. Fine tuning of noise magnitudes. Once noise compo- 
nents have been identified in frequency and magnitude a fine 
tuning estimate is obtained by scanning the modulation fre- 
quency and fitting to the full spectrum to Eq. (J3J) (red line). 
Here we use n = 11 equidistant pulse sequence. 

What modulation suits best for spectrum estimation? 
Ideally, one would require a Fourier-limited sine around 
the modulation frequency, as used in [21]. It is, how- 
ever, more convenient experimentally to use a stream of 
7r pulses. For the purpose of noise spectroscopy Yuge 
et. al. [3T] suggested an equidistant pulse scheme while 
Cywihski et. al. [TS] suggested the Uhrig [5] scheme. Pre- 
viously, these schemes have been compared experimen- 
tally with trapped ions |22j and solid state NMR [34 in 
the context of optimized dynamical decoupling. Here we 
use both schemes for the purpose of spectral estimation. 
A comparison of typical modulation spectra of the two 
is shown in Fig. [4^i and [4]b; in both cases 19 pulses are 
used for a total experiment time of T — 66.7 ms. 

In terms of spectral filters, the equidistant case is 
more intuitive since its peak is centered around f mo d = 
(n + 1)/2T while in the Uhrig case, the spectral peak 
is shifted towards lower frequencies. Moreover its spec- 
tral content at high frequencies is much greater than the 
equidistant case. However at low frequencies, Uhrig mod- 
ulation decays faster. It appears that Uhrig modulation 
should be more suited when one wants to measure dis- 
crete high frequency components in the presence of un- 
wanted low-frequency noise. Otherwise equidistant mod- 



TABLE I. Fine tuning of field magnitudes using equidistant 
and Uhrig modulation. 
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ulation is simpler to use and interpret. 

The intcrprctational simplicity of equidistant modula- 
tion is pronounced in the two-dimensional scan shown 
in Fig. Hp. For each fixed experiment duration, T, a 
phase scan is shown as a column in figure [4p. The con- 
trast of each column is obtained from a fit procedure and 
displayed as a single data point in Fig. [4p correspond- 
ingly. The noise spectral peaks can be easily identified 
at 150 Hz, 200 Hz and 250 Hz. The spectral shape of 
the Uhrig modulation (Fig. |4Jd) mixes nearby frequency 
components of the noise, as seen in figure [4}i and [4]:. 

Quantitatively, with our noise profile, both modulation 
schemes performed equally well and show agreement in 
the extracted spectrum (figure |4^,f). This is summarized 
in the table |l| The Uhrig scheme reveals information on 
an additional frequency component, Biqqhz, due to the 
shift of its spectral peak to lower frequency noise as com- 
pared with the equidistant case. 

In conclusion, we have developed and used an an- 
alytical model for spectral noise estimation of non- 
perturbative, non-Gaussian, discrete dephasing noise. 
We used our scheme to extract the magnetic field noise 
spectral components at the power line harmonics in var- 
ious real lab scenarios. In fact, this enabled us to cali- 
brate our magnetic field compensation system to reduce 
noise components to the fiG level, thereby reaching 1.4 s 
of coherence with a Zeeman sensitive qubit 30J . We 
expect this model to be useful for other discrete noise 
scenarios. One example is spontaneous a-oscillations in 
brain activity measured using magnetoencephalography 
(MEG) [35] ■ Another example is the study of decoher- 
ence of a single NV center induced by a finite number 
of 13 C nuclear spins in diamond 36J or the discrete me- 
chanical resonances of a cantilever coupled to an NV cen- 
ter [31] ■ 
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FIG. 4. Comparison of two modulation schemes: equidistant and Uhrig. (a) The spectrum of equidistant modulation centered 
at 150 Hz with 19 pulses, (b) Same as a, for the 20 pulses Uhrig scheme with the same total duration of 66.7 ms. The 
modulation spectrum has much more energy in high frequencies, and its peak is shifted to 109.6772. (c) Each column is a scan 
of 4> r j at n = 19 equidistant pulses and a fixed T. The total experiment time T is varied from column to column and the 
corresponding modulation frequency is shown in e. (d) Same as c, for an Uhrig scheme of 20 pulses, (e) Each point corresponds 
to the fringe contrast of corresponding column in c. Clear peaks are pronounced in 150 Hz, 200 Hz and 250 Hz. Red line is a 
fit to Eq. |3|. (f) Same as e, for the Uhrig scheme. As the modulation mixes neighboring frequencies, it is harder to identify 
spectral peaks compared to the equidistant case. The fit indicated by the red line yields field magnitudes consistent with the 
equidistant modulation fit in e. 
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